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ABSTRACT 

The expl oration of mean-field gal actic dyna mos affec ted by a galactic spiral pattern, 
begun in IChamandv et al.1 (|2013l hereafter iPaper J ) with numerical simulations, is 
continued here with an asymptotic solution. The effect of the spiral pattern on the 
dynamo considered is the enhancement of the a effect in spiral arms (which may 
overlap the gaseous spiral arms or be located in the interarm regions). The axisym- 
mctric and enslaved non-axisymmetric modes of the mean magnetic field are studied 
semi-analytically to clarify and strengthen the numerical results. The mean-field dy- 
namo model used generalizes the standard theory to include the delayed response of 
the mean electromotive force to variations of the mean magnetic field and turbulence 
parameters (the temporal non- locality, or r effect). Good qualitative agreeme nt is ob- 
tained between the asymptotic solution presented and numerical solutions of iPaperll 
for a global rigidly rotating material spiral (density wave). We find magnetic arms 
displaced from the a-arms away from the co-rotation radius, so that the ridges of 
magnetic field strength are more tightly wound than the a-arms. Moreover, the effect 
of a finite dynamo relaxation time r is to phase-shift the magnetic arms in the di- 
rection opposite to the galactic rotation. This mechanism can be used to explain the 
observed phase shifts between magnetic and material arms in some galaxies. 

Key words: magnetic fields - MHD - galaxies: magnetic fields - galaxies: spiral - 
galaxies: structure - galaxies: ISM 



1 INTRODUCTION 

■ Disc galaxies typically contain regular (or large-scale or 
mean) magnetic fields of strength ~ 1 — 10 /iG. In many 
cases, the strength of the non-axisymmetric components 
can b e comparable to that of the axisymmetric compo- 
nent l|Fletcherl liuToh . The non-axisymmetric part of the 
field takes the form of ' magnetic spiral arms', wherein the 
regular fiel d is enhanced l|Beck et al.lll996l : IShukurovll2005l ; 
lBeckll2012f ). Magnetic arms cannot be easily accounted for 
by mean-field dynamo theor y if the underlying disc is as- 
sume d to be axisymmetric (Ch amandy et all 1201 3l . here- 
after jPa^iO) . Moreover, they seem to be related (though 
in a non-trivial way) to the material (gaseous) spiral arms. 
Hence the need to develop a theory which relates the non- 
axisymmetry of the regular magnetic field to that of the 
underlying disc. This was the general aim of lPaperll . where 
we approached the problem from a numerical standpoint and 
considered not only the linear growth of the field (kinematic 
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regime) but also the non-linear saturation phase. Here we 
take an analytical approach and therefore restrict the inves- 
tigation to the kinematic regime. We are motivated, in part, 
by insig hts that were provided by analy tical tre atments in 
the past (|Ruzmaikin et a l. 1988, he reafterlRSS88l, for the ax- 
isymm etric case and lMestel fc Subramanianl f99fl . hereafter 
MS9f, for the non-axisymmetric case). The latter work was 
focused on how the |m| = f (bisymmetric) azimuthal modes 
can be generated in a galactic disc with a two-arm spiral pat- 
tern. Nu merical models (e.g. [Moss f 99§j; iRohde et al.lll999l : 
IPaper if ) show, however, that m = (axisymmetric) and 
m | = 2 (quadrisymmetric) modes tend to be more preva- 
lent in such discs (or m = and \m\ — n for discs with 
n spiral arms). An analytical model that can explain such 
modes is therefore needed. 

Our aim in this paper is thus a systematic interpretation 
of th e key numerical results of IPaper il (| Chamandy et all 
l2013h by solving a suitably simplified but essentially the 
same mean-field dynamo equations with an approximate 
semi-analytical method. Here we consider both axially sym- 
metric and enslaved non-axisymmetric modes of a kinematic 
dynamo. The enslaved modes are those that have the same 
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growth rate as the leading axisymmetric mode; they occur 
because of deviations of the galactic disc from axial symme- 
try, e.g., due to a spiral pattern in the interstellar gas. For 
a two-armed spiral, the enslaved modes include the \m\ = 2 
modes that corotate with the spiral pattern, as well as other, 
weaker, even-m corotating modes. We leave the \m\ = 1 
modes and other odd-m corotating modes for a fut ure wor k. 
Our model differs from earlier analytical study of lMS9ll in 
that here we: 

(i) provide an asymptotic treatment of the m = and 
even-m modes forced by a rigidly rotating two-arm spiral; 

(ii) incorporate the minimal-r approximation (MTA) and 
explore the effects of a finite dynamo relaxation time. 

The plan of the paper is as follows. The basic equations 
and mathematical approach are summarized in Section [2] 
In Section [3] we present an asymptotic solution for the ax- 
isymmetric and enslaved non-axisymmetric magnetic modes 
which inhabit a disc containing a global, steady, rigidly ro- 
tating spiral density wave. We then compare our results with 
numeric al solutio ns in Section [3] Conclusions and discussion 
for both iPaper 11 and this paper can be found in Sections [5] 
and [6] 



namical equation for £ : 

§ = I(cB- w VxB)-f. (2) 

In the kinematic limit, and for isotropic and homogeneous 
turbulence, the turbulent transport coefficients are given by 

a = — \t c "u ■ V X u, r/ t = \t c u 2 , (3) 

with r c the correlation time of the turbulence and r the 
relaxa tion time, found in numerical simula ti ons to be close 
to r c (jBrandenburg fc Subramanianll2005al fbl 120071) . We re- 
tain c T — t /t c as a, free parameter, but set it to unity in 
illustrative examples. For the sake of simplicity, we take r 
to be independent of scale and constant in space and time. 
Eq. ((2]), when combined with the mean- field induction equa- 
tion, leads to Eq. (fjQ). 

The more standard mean-field dynamo equation 

?Tf5 

— = V x ([/ x B + o.B) + ?jt V 2 B, (4) 

is recovered from Eq. Jl} in the limit r — > with c T = 1. 
Below, we refer to the approach leading to Eq. Q as the 
'standard theory'. 



2 THE MEAN-FIELD DYNAMO AND THE 
TELEGRAPH EQUATION 

The basic equation solved here is a generalization of the 
standard mean-field dynamo equation, and is derived in 
IPaper J , 



d 2 B OB 

T dt 2 + dt 



(1) 



+ V x (U x B + c T aB) + c T r/ t V 2 B. 

Here we use the same (standard) notation as in IPaper J . 
We have neglected the terms rr]X7 2 (dB/dt) and 77V 2 B on 
the right-hand side because of the high Ohmic conductivity 
of the interstellar gas, 77 <C rjt- Following the usual mean- 
field approach, the velocity and magnetic fields, U and B, 
have each been written as the sum of an average part and a 
random part, 

B = B + b and U = U + u, 

where overbar represents the ensemble average but for prac- 
tical purposes spatial averaging over scales larger than the 
turbulent scale but smaller than the system siz e can be used 
(|Germanolll992l ; lEvinkll2012l ; iGent et alj|2012h . 

Equation fTJ arises when the response time r 
of the mean electromotive force (emf) £ — u x b 
to changes in the mean fi eld or small-scale turbu- 
lence is not negligible (e.g. iRheinhardt fe Brandenburg! 
l2012h . The t effect can have various interesting astro- 
phys ical implications (see, e.g., iHubbard fc Brandenburg 
120091 . which also contains a brief review of applications 
from the literature). A convenient approach to allow 
for this effect is the minima l- r approximation ( MTA) 



llRogachevskii fc Kleeorinl |200C|; iBlackman fc Field! |200 



3 ASYMPTOTIC SOLUTIONS FOR 
GALACTIC DYNAMOS 

Many disc galaxies have a spiral structure which causes de- 
viations from axial symmetry in both turbulent and regu- 
lar gas flows, l eading to non-axisymmetric magnetic fields. 
Here we follow |MS9lJ by assuming that the a effect is modu- 
lated by the spiral pat tern. The nature of such a modulation 
is lar gely unexplored l|Shukurovlll99i ; I Shukurov fc Sokolofj 
1998), but it is reasonable to assume that a is enhanced 
along a spiral perhaps overlapping with the gas spirals. The 
spiral patterns observed in the light of young stars are be- 
lieved to represent a global density wave rotating at a fixed 
angular frequency, or possibly a superposition of such waves, 
with diff ering azimuthal symmetries and angular frequen - 
cies fe.g. IComparetta fc Quillenll2012l : iRoskar et al.l 12012'). 



(e.g. Dobbs et al. 


l20ld; Sellwood 201ll;lQuillen et al.ll201ll; 


IWada et all 12011 


: iGrand et al.l 120121) whose effect on the 



we consider an enduring spiral with n arms and a constant 
pattern speed f2. 



3.1 Basic equations 

Since we are interested in non-axisymmetric magnetic fields 
that corotate with the spiral pattern, it is preferable to 
work in the corotating frame w here th e dynamo forcing is 
independent of time. Following MS91, we carry out a co- 
ordinate transformation from the inertial cylindrical frame 
E = (r, cj>, z, t) (with disc rotation axis as the z-axis) to the 
frame E = (r, <j>, z, t) rotating with the pattern angular ve- 
locity Q (assumed constant in both time and position): 



iBrandenburg fc Subramanian 2005a ). which leads to a dy 



= (j> ~ Qt, r = r, z — z, t = t, 
B r = B r , £>0 = B^,, B z = B z . 



(5) 
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so that 



d 



The coordinate transformation shifts the mean velocity, 
U$ = 17 ^ — rfi, but the random velocity is left unchanged, 
u = u, as are both mean and random magnetic field. There- 
fore, 5 = a and 5/t = rjt ■ 

We use the thin-disc approximation, d/dz 2> <9/9r 
and, eventually, consider tightly wound magnetic spirals, 
d/dr d/rd(j>- Thus, we consider magnetic fields whose ra- 
dial scale is asymptotically intermediate between the scale 
height of the galactic disc (of order 0.5 kpc) and its radial 
scale length (of order 10 kpc) . We also adopt the atj-dynamo 
approximation where the rate of production of the azimuthal 
magnetic field by the a effect is negligible in comparison with 
the effect of the differential rotation. The mean velocity field 
is taken to be axisymmetric, purely azimuthal, and constant 
in time, 



U — roj(r)<j>, 



(7) 



where </> is the unit azimuthal vector. In the rotating frame, 
the gas angular velocity is w(r) = U}(r) — £1. The magnitude 
of the rotational velocity shear is quantified with G(r) — 
rdto/dr — rdu/dr. Applying Eqs. J5]) and to the r- and 
(^-components of Eq. {TJ, and dropping tildes for ease of 
notation (except for on ui), leads to the following equations 
in the corotating frame: 



CB r 



9 I n 



oz 

+ c T r/t ( V 2 B r + 



1 d 2 B r 
r 2 dij> 2 



2 dB 



1+T 



dt 



nr-) {GB r 



+ c T r, t V B^ + 



d<p 



2 dB, 



(8) 



(9) 



where 



and 



dt 



fir 







d<p J \dt 



d 



+ 0J- 



d 2 X d 



+ 



dr 



ld_ 

r dr 



(rX) 



These equations agree with Eqs. (2.2) and (2.3) of lMS9ll in 
the limit r — ¥ with c T = 1. Since, in the thin-disc approxi- 
mation, both equations do not include B z , the equation for 
B z can be replaced by the solenoidality condition, 



1 d - 1 dBj, 8B Z 
--^-(rBr) H -^j- + 

r or r oq> oz 



= 0. 



(10) 



We emphasize that the governing equations acquire 
additional terms when transformed to the rotating frame 
whenever r is finite. On the contrary, the standard mean- 
field dynamo equation @ is invariant under the transforma- 
tion to a rotating frame. This can be understood physically 
in the following way. In the standard theory, £ at time to 
and position x depends on B and the turbulence param- 
eters at the same time and position. This is not true in 
the MTA, where £(to,x) depends on the history of B(t,x), 



a and rjt over, roughly, to — t < t ^ to- [Here 5 and rj t 
refer to the more general version of Eq. ([2]), wherein we 
have made the standard assumption 5 w a/r c , rj t ~ Tjt/r c 
l|Brandenburg fc Subramanianll2005al ).] Thus, £ (to,x) is af- 
fected by how fast variations in the field and (a, %) are 
sweeping past the position x. In other words, the dynamo 
equations are no longer invariant under the transformation 
to a rotating frame. 

3.2 Approximate solution 

When the coefficients of Eqs. (J5J and © depend on the az- 
imuthal angle (j> as exp(im0) with certain m, their solutions 
can be represented in the form 



m — — oo 



a m (r, z) 



b m (r, z) 

so that Eqs. ((8]) and ((9} reduce to 



exp (im<t> + Tt), (11) 



E 



m— — oo 



Ca m + c T -^-(abr, 

Oz 



-c T r)t (V a m — m — 2im--^ 



Cb m — (1 + Ft — imQ.T)Ga r , 



c T r]t V b 



2 bm , ; , . Clm 

m — + 2im— - 



0, (12) 



= 0, (13) 



where 



C = (1 + Tt - imf2r)(r + imw). 



We note that a is, in general, a function of r, <f> and z, 
whereas cj and G are assumed to be functions of r alone. 
When obtaining Eq. (|13|) from Eq. ([9|, the summation can 
be dropped because the only (^-dependence that occurs in 
this equation is in the common factor exp (imcj)), whereas a 
depends on <f> in Eq. (|12|) . 

To make further progress, we model a as an n-armed, 
rigidly rotating global spiral (thus, in the rotating frame it 
is time- independent), 



= ao[l + e Q cos(n(f> — nr)] = ao+a, 



in4> 



+a- 



where 



(14) 



(15) 



where ao may depend on r and z, and asterisk denotes 
complex conjugate. With this convention, n < describes 
a trailing Archimedean spiral. We choose an Archimedean 
rather than a logarithmic spiral in order to simplify the cal- 
culations; this choice does not affect our conclusions. The 
requirement that each coefficient of Eq. ()12[) vanishes gives 

r ( ^2 2 dm n . b m 
La m —c T r) t V a m — m — 2im— 



~ —c T — (aob m + a n b m - n + a_ n b m+n ). 
oz 

This implies that modes with m = — In, — n, 0, n, 2n, ... 
are coupled to each other, and likewise modes with m = 
— 2n + 1, —n + 1, 1, n + 1, 2n + 1, etc. However, these 
sets of modes are decoupled from one another. Modes that 
are coupled to the generally dominant m — mode grow 
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along with it, and are thus called 'enslaved' modes. One 
naturally expects the lowest-order enslaved modes of order 
n to dominate the higher-order enslaved modes, as the for- 
mer are forced dir ectly. Th is was in fact borne out in the 
numerical work of IPaperll Therefore, the m = ±n modes 
will be the dominant non-axisymmetric modes for a galaxy 

with an n-armed spiral pattern. 

Taking n = 2 ( a two- armed a- spiral), lMS9ll and 
ISubramanian fc Mestel l| 19931 . hereafter ISM93h truncated 
the series of coupled equations by neglecting \m\ ^ 3. Thus, 
only the m = ±1 are coupled and likewise m = is cou- 
pled to m = ±2. Here we tru ncate the series in the same 
way; the numerical solutions of |Paper"H confirm that this is 
an excellent approximation for much of the parameter space 
(however, see below) . The case of m = ±1 has already been 
considered in MS91. Thus, we focus on the m = 0, ±2 case 
here, for which the governing equations are 

~2 d 

(1 + r-r)rao — c T ?7tV a = —c T — (aobo + «2b-2 + a-262), 



(1 + rv)rb - crJjtV 2 io = (1 + rv)Ga , 

while, for m = ±2 and n = 2, 

f~2 . «±2 ^_ . .&±2 

La±2 — c T r/ t I V a±2 — 4— — =p 4i — j- 

= — c t -^-(qo&±2 + a±2bo), 
oz 

= (l + rr=F2ifir)Ga ±2 , 



(16) 
(17) 



(18) 



(19) 



where 



equations are approximated as algebraic equations, works 
well for the case where the vertical turbulent diffusion is 
much stronger than the radial turbulent diffusion, i.e., for 
solutions of a large radial scale (wavelength) in a thin disc. 

In the no-z approximation, the vertical diffusion terms 
are approximated as 



d 2 a±2 
dz 2 



a±2 
Ah 2 ' 



d 2 b ±2 
dz 2 



Ah 2 ' 



where h is the disc half-thickness. Further, the terms con- 
taining a can be approximated as 

g 2 

— (a b + a 2 b-2 + a.-ibi) — —r(aobo + a 2 b- 2 + a-262), 
dz Tih 

-^-{a>ob±2 + a±2bo) ~ — (ao&±2 +a±2&o)- 
oz Tin 

For the radial diffusion terms, we start by representing 
a±2 and b±2 as 



a±2 



i(r)e Ti( ' W , 



b±2 = b(r)e 



(20) 



where the minus sign in front of a is introduced for future 
convenience, and where we have removed the z-dependence 
of a and b due to the fact that we are now working in the 

no-2 approximation. Note that \a%\ = | a,_ 2 1 and |&2| = | & 2 1 , 

whereas the phases of 0,2 and a_2 (or 62 and b-2) must 
be equal in magnitude and opposite in sign, since B r and 
B$ are real. We then apply the WKBJ-type tight-winding 
approximation, assuming that the amplitude varies much 
less rapidly with radius than the phase, i.e., 



da 








db 




-d(h 




< 


a— 






< 




dr 


dr 




dr 


dr 



The radial diffusion terms are then approximated by 
d 



L = (1 + Tt =F 2inr)(T ± 2vJj). 



dr 



Id 

--7-( m ±2 

r dr 



{9'a ± iO'a) a± 2 , 



(21) 



3.2.1 Non-axisymmetric modes under the no-z and 
WKBJ approximations 

We note that for the even modes m = 0, ± 2 there are six 
coupled partial differential equations. lMS9ll solved the cor- 
responding simpler problem for the odd (m = ±1) modes, 
which involves four coupled PDEs, using WKBJ methods, 
but they did not do the same for the even modes. Thus 
we would like to get analytical insights into these perhaps 
more important modes. Therefore we proceed as follows: 
Our plan is to first solve for a±2 and b±2. We do this by 
using the no-z approxim ation for the z-derivatives {SM93; 
lMossill995l ; |Phillips|[200l1) , and a WKBJ-type approximation 
to handle radial diffusion. This yields algebraic equations 
which can be solved to obtain explicit expressions for a±2 
and b±2. We then substitute them into the m = equa- 
tions, (|16p and (I17[l . and apply the no-z approximation to 
obtain Schrodinger-type differential equations in r with the 
effective potential V(r). Bound states in this potential cor- 
respond to exponentially growing solutions, and these are 
obtained using the WKBJ approximation. This is done it- 
eratively using numerical methods, so that the method is 
semi-analytical rather than analytical. Thus we avoid solv- 
ing the six coupled differential equations (|16|l - (|19|) . We shall 
see that this procedure, in which four of the six differential 



d 
dr 



1 d 
r dr 



(rb±2 



t±i%)b ± 2, 



(22) 



where prime stands for d/dr. It should be noted that in 
the region where the relative strength of non-axisymmetric 
modes with respect to axisymmetric modes is expected 
to be the largest (near the co-rotation radius), we expect 
1/r <C I d/dr \ since the co- rotation radius is typically com- 
parable to the size of the galaxy. Thus, we have dropped 
terms proportional to 1/r. Since turbulent diffusion terms 
containing the (^-derivatives are proportional to 1/r 2 , they, 
too, can be neglected. 

As shown in Appendix [XJ solutions of Eqs. (|18[) and 
H19|) can be represented as 



a±2 = 



\Do\b e a A 2 + B 



exp[=Fi(Kr + /3 + /?)], (23) 



b±2 = 



\D \boe a 



2 VA 2 + B 



A 2 + B 2 
1 

exp [=pi(«r + f})\, 



(24) 



where A, B, A and B, as well as (5 and j3, are functions of 
position defined in Appendix [X] Here the dynamo number 
is defined as 



Do = 



a Gh 3 _ a Gt 2 d 



<0, 



(25) 
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where td = h 2 /rjt is the vertical turbulent diffusion time 
scale. 

Using Eqs. (|20|) . we now make the identifications 



fir + p + P, 



% = nr + P, 



and so 



k+p', e': = p"+p", 



--ft'. 
(26) 



Clearly, P and P depend on A, B, A and B (Eqs 
and lA8|) , which, in turn, depend on the derivatives of P and 
P through 9' a , 9' b , 9'a and 6£ (Eqs. |H H3 \M and IATO)) . 
This means that, unless P' and /?' are completely neglected, 
it is not possible to obtain an analytical solution. We will 
see however that, if the terms involving P' , P' , P" and p" 
are small compared with the turbulent diffusion terms that 
do not vary with r, then P and P can be determined semi- 
analytically, by successive iteration. 

The need for the iterations can be understood as fol- 
lows. Firstly, if P' and P' were zero, then the radial phase 
would be «r, which is the same as the radial phase of the a- 
spiral. Thus, the magnetic arms would be as tightly wound 
as the a-arms. In fact, the extra contributions to the radial 
phase P + P for B r and P for B,/,, are finite, and turn out 
to have the same sign as k. This implies that the magnetic 
arms (defined as the set of positions where B is maximum), 
are more tightly wound than the material arms that pro- 
duce them. This is because the enhancement of field due to 
the enhanced dynamo action within the a-arm is advected 
downstream with the differentially rotating gas. Ultimately 
a balance is set up between the competing effects of prefer- 
ential magnetic field generation in the a-arms and advection 
downstream, resulting in stationary magnetic arms that are 
more tightly wound than the a-arms, and cross them near 
the co-rotation circle. Moreover, this balance, which deter- 
mines the pitch angle of the magnetic arms (ridges), and thus 
the values of P and P, involves many quantities, including 
the derivatives of P and P through the radial diffusion (as 
seen in the equations of AppendixfSJ. This self-dependence 
results from the fact that the radial diffusion is stronger 
when the magnetic arms (ridges) are more tightly wound, 
as can be seen from Eqs. (|21[) and (|22[) . Thus, the radial 
diffusion feeds back onto itself, forcing us to iterate to obtain 
the functions P and p. 



using Eqs. (127[) . we can write the overall solution as 



B r 



-aq 



A 2 + B 1 
A 2 +B- 



e a \D \cos(2<f> - nr ~ P - P) 



B$ =bq 



1 



=e a |-D | cos(20 - nr - P) 



(28) 
(29) 



VA 2 + B 2 

Magnetic field lines of the mean field are expected to be trail- 
ing spirals since G < 0, which implies b/a < 0. Therefore, 
the overall phase difference between B r and B^ has magni- 
tude P [without the minus sign introduced in front of a in 
Eq. I|20p. there would have been an extra phase difference 
7r]. The pitch angle of the magnetic field, defined as 



Pb = arctan 



B r 
B„ 



(30) 



is independent of q. 

Each component of the mean magnetic field consists of 
an axisymmetric part and that having the m = 2 symmetry, 

Bi = Bf ] +B? ) 008(2(0 - &)], 

where i = r,(j>. The magnitude of the r-component of the 
field is maximum where 20 — nr — P — p = 0, while the 
magnitude of the 0-component is maximum for 2(f> — nr — 
P — 0. On the other hand, the magnitude of a is maximum 
where 2(f) — nr — 0. Therefore, the (7r-fold degenerate) phase 
differences between the r and <f> components of the magnetic 
field and the a-spiral are given, respectively, by 



A r = 



P + P 



2 ' 



3.2.3 Growth rates of non- axisymmetric magnetic modes 

We are now prepared to obtain the growth rate V. For this 
purpose we substitute Eq. (|27p into the no-z versions of 
Eqs. I)16p and (JTTj) , and, in Eq. (|16[) . we also substitute 
Eq. (|A1[) for b±2. In addition, we make use of the relation 

2 2/. 

020-2 = a e a /4, 

obtained from Eq. (I15[) . Equations (|16[) and (|17[) then reduce 
to 

>g)' V 



(1 + Tr)r - c T 



2c T b 
— rtto - 
Tin a 



1 + \£2lA {X2 + x _ 2: 



c T r]t 



4h 2 ' 



3.2.2 The global solution 

Once a±2 and b±2 are obtain ed, we c an solve E qs. (|16|l and 
(TIT)) for oo and b . Following IRSS88I and lMS9ll by invoking 
the thinness of the disc, we factorise the solution into the 
local (5 and b) and global (q) parts, 



a (r) = aq(r), b (r) = bq(r). 



(27) 



Generally, the local solution depends on z, i.e., a(z) and 
b(z). However, the ^-dependence has been removed using 
the no- ,z approximation, and the local solution is a vectorial 
constant. Substituting Eqs. J23} and J24} into Eq. (fTTj) . and 



(i+iv)r- 



(rq)' 



(1 + Tt)G= - cvijt 



Ah 2 ' 



where X±2 are defined in Eq. (|A2[) . Note that e a — cor- 
responds to an axisymmetric forcing of the dynamo. It is 
convenient to introduce a new quantity 7(r) defined as 



7 (r) = (l + rr)r-c T ^ 



(rq)' 



whose significance will shortly become evident. Now the 
equations can be separated into a global eigen-equation for 
q(r), which only contains derivatives with respect to r, and 
two local equations for a and b, which generally only have 
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derivatives with respect to z but become algebraic equations 
under the no- z approximation. 

The global equation is given by 



(rq)' 



(i + rv)r g = 79 + c T Tjt 

while the local equations can be rewritten as 

\D \el 



(31) 



7 + 



2a ~ 
c T — —b 
nn 



1 



(X2+X-2] 



=(.1 -LVKV/i. 



(32) 
(33) 



It is now clear that 7 is the local (at a given r) growth rate 
of the dynamo modes in the slab approximation. 

These homogeneous equations in a and b are easily 
solved to yield b/a and 7. From Eq. (|33|1 . we obtain 

6 = (l + Vr)Gt d 

Then the vanishing of the determinant of the coefficients of 
Eqs. (32) and J33) yields 



?SL(1 + Vt)\D \ 

7T 



1 + 



e%\D \A 



2(A 2 + B 2 



(34) 

where the positive sign in front of the square root provides 
growing solutions, 7 > 0. 

With e a = (axisymmetric forcing) and 7 = (neutral 
stability of the mean magnetic field), we obtain an estimate 
for the local critical dynamo number, i.e., the minimum mag- 
nitude of the dynamo number for the field to grow rather 
than decay, 



-Do.cr = 



9.6c T 



32(1 + IY) 1 + IY 

It is slightly larger than the Do.cr = 8, obtained by the 
numerical solution of the z-depende nt versio ns of Eqs. ()32[) 
and (EH) for c T = 1, r = e a = (RSS88); we note that 
Vt < 1. 

For e a = 0, the magnetic pitch angle, given by Eq. (|30[) . 
can be expressed as 



tanps.o 



2c T 



Tv(l + Fr)h\G\' 



(35) 



A similar result can be obtained from perturbation theory 
(|Sur et alj|2007h . With F. Krause's law, 

ato = l 2 uj/h, (36) 

and G = ~~u) (flat rotation curve), this becomes 



I / 2c r 

tanpBi0 = -_y ___ . 

The azimuthal average of the pitch an gle in the kinematic 
regime obtained numerically in lPaper H where Eq. ((36j) was 
used, agrees very well with this formula. This can be seen by 
using the relevant parameters from that paper in the above 
equation, and then comparing with Fig. 8d of that paper. 



However, here we shall take Qo and h constant with radius, 
so that we have 



tanp fl ,o 



2c T 



on 



tt(1 + Vt) huj' 



(37) 



so that \pb,o\ now increases with radius, as uj is a decreasing 
function of r. 

Now, substituting Eq. (|34p into Eq. (|3ip . we obtain a 
Schrodinger-type equation for q{r), 

-\7 2 r q + V(r)q(r) = Eq(r), (38) 

where V? = (rq )' jr is the radial part of the Laplacian in 
cylindrical coordinates, with prime denoting d/dr. Here, the 
'potential' is given by 



V ■■ 



r 2 Crh 2 



and the 'energy' eigenvalue, by 



1 + 



el\D \A 



2(A 2 + B 2 ) 



E : 



(1 + rr)rt d 

Orh 2 



(39) 



(40) 



Note that bound states in the 'potential' V(r) with E < 
correspond to growing modes with V > 0. We solve this 
equation using the WKBJ theory: in terms of a scaled vari- 
able x, introduced via r = he x (with h = const, — 00 < x < 
00 and dx = dr/r), Eq. (|38|) reduces to 

d 2 q 

— 2 +P(x)q{x) = 0, 



dx 2 



where 



he [E ■ 



V(x)}. 



The boundary conditions are q — at r = and r — > 00 (or 
x — > ±00). Suppose that p(x) = for x = x± (corresponding 
to r±), with X- < x+ (for our purposes, it is sufficient to 
consider the case of just two roots). The standard WKBJ 
theory yields the quantization condition 



-dr = 



~[E-V(r)] 1/2 dr = ^(2k + l), (41) 



with k = 0, 1,2, . . which can be used to obtain the growth 
rate V. Note, however, that V enters both the 'potential' and 
the 'energy', so the WKBJ approach has to be supplemented 
with an iteration procedure to converge on the correct value 
of r, in addition to the iterations over /3 and /3 discussed 
above. 

As we will see below, the ^-dependent term in the a 
effect leads to a deepening of the potential well near the 
co-rotation radius, allowing non-axisymmetric modes that 
corotate with the a-spiral to exist there. Outside this re- 
gion (i.e., far from the co-rotation radius) axisymmetric 
modes will dominate, while non-axisymmetric modes will be 
weaker. Non-corotating non-axisymmetric modes (which ro- 
tate at approximately the local angular velocity at the radius 
where their eigenfunctions are maximum) are also exp ected 
to exist, but to be sub-dominant everywhere (|MS91). The 
asymptotic solutions developed here apply to the corotating 
magnetic mod es in the kinematic regime, but the numerical 
simulations of |Paper"H are, of course, not restrictive. 



© 2012 RAS, MNRAS 000,[ljfl4] 



Galactic spiral patterns and dynamo action II 7 




Figure 1. (a) The potential (130 P of J38D for r = 21 /u (thick, black, solid curve) and t — > (thin, red, solid curve), near the co-rotation 
radius (dashed vertical line). The 'energy' eigenvalues l|40|l of l|38|l are represented by horizontal lines inside the wells. The axisymmctric 
part of the potential is also shown as a dashed line in both cases, (b) The solid curves represent the ratio of non-axisymmetric to 
axisymmetric parts of B^ at the azimuth cj> = <f> cot = /tr cor /2, while the dashed curves show the corresponding magnitudes of the ratios 

±B ( ^ ) /B ( ° } . The dash-dotted curve shows the variation of a(r, <fr C or)- (c) The phase differences, A0 (solid) and A,- (dashed), between the 

components of the magnetic spiral arms and the o-spiral arms. The radius r max where the envelope of B? /B? is maximum is shown 
as a dotted vertical line of the appropriate colour and thickness, (d) Pitch angle pg, shown with the opposite sign, at cj> = (j> cor = Kr cor /2 
(solid). Also shown is — pg for the purely axisymmetric case (dashed). 



3.3 Illustrative example 

We determine V and the properties of the mean magnetic 
field by varying V iteratively until the quantization condi- 
tion (|41[) is satisfied for k — (the fastest growing mode). 
In addition, there is the complication that /3, j3, and hence 
9 a , Ob, A, B, A, B, etc., are not a priori known (see Sec- 
tion l3.2.T1 and Appendix[X]| ■ As mentioned above, we resolve 
this by iterating these variables as well. That is, we start 
with P'(r), fi'(r) and their derivatives equal to zero, and at 
each iteration, we obtain these functions from Eqs. (|A5|) and 
(|A8|I to use them as input for the next iteration. The iter- 
ations converge if the terms involving /?' and ft' are small 
compared with the turbulent diffusion terms that do not 
vary with r. In the example below, /3' and /3' are compara- 
ble to k near the co-rotation. However, it can be seen from 
Eqs. (|A6p. (|A7[1 . (|A9f) and (|A10p that the terms involving 
the radial derivatives are not very important if the radial 
diffusion is much weaker than the vertical diffusion, that is 
if 8' a 2 , 9' 2 <C ty 2 /4h 2 , which may permit the iterations to con- 
verge. Therefore, in the example below, we choose the disc 
half-thickness h to be small enough (somewhat smaller than 
the standard value) that the iterations do in fact converge. 

In order to determine the potential and growth rates 
of the corotating non-axisymmetric modes, we must first 
specify the functional forms of the scale height h(r), an- 
gular velocity oj(r) and ao(r). For simplicity, we choose 



h — 0.1 kpc = const. We also take % = ul/3, where u 
and I are the rms turbulent velocity and scale, and adopt 
l = h/A = 0.025 kpc, and u = 12/if." 1 = 5kms~\ These val- 
ues are different from those usually adopted, u — lOkms -1 
and I = 0.1 kpc. Our choice is motivated by the desire to 
improve the convergence of the iterations since our aim here 
is mainly to illustrate the solution procedure and to discuss 
the qualitative properties of the eigen-solutions. However, 
we take care to keep the key parameter, the vertical turbu- 
lent diffusion time, at id = 0.24 Gyr, similar to (about twice 
smaller th an) the standard value. 

As in lPaperll . we use Brandt's rotation curve, 

wo 

U ~ [l + (r/ rul ) 2 } 1/2 ' 

with = 20h = 2 kpc, and ujq set to yield the circular 
rotation speed U $ = 600/it^ = 250kms _1 at r = lOO/i = 
10 kpc. We choose the co-rotation radius to be located at 
r CO r = 100ft = 10 kpc, which corresponds to a pattern speed 
n = 6^ = 25kms- 1 kpc _1 . 

Furthermore, we take «o = 2.4/itJ 1 = lkms" 1 at all 
radii. We also take n = — O.l/i. -1 = — lkpc -1 and e Q = 1, 
the latter chosen to be large to maximize the strength of the 
non-axisymmetric modes. We also adopt cv = 1. 

The value of r is estimated as r w 21 /u ~ 2/ 2 /3r;t ~ 
2l 2 t d /3h 2 = td/24 = 9.8 Myr. We also consider solutions for 
t — > 0, both for reference and also because r could be much 
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smaller than our estimate of 21 /u in some galaxies (though 
it could be much larger in others). 

When describing the geometry of the magnetic modes, 
it is useful to have in mind that magnetic lines have 
the shape of a spiral with the pitch angle defined in 
Eq. (|30p . These spirals have to be distinguished from the 
spirals along which magnetic field strength is maximum. 
These can be called magnetic arms or magnetic ridges , 
and have a different pitch a ngle jBarvshnikova et al.lll987l ; 
iKrasheninnikova et al.lll989h (see also Fig. VII.8 of lRSS88l) . 

The converged results are shown in Fig. [1] for r — > 
and r = id/24 = 9.8 Myr. (Solutions for intermediate val- 
ues of r are intermediate between the two solutions il- 
lustrated.) Panel a shows the potential V(r) of Eq. (139[) . 
along with the 'energy' — (l + 'jT)'yta/h 2 of the fastest grow- 
ing eigenmode. Importantly, rue to the non-axisymmetric 
part of a, the potential has a local minimum near the 
co-rotation radius. Therefore, corotating non-axisymmetric 
modes can be preferentially excited there. The finite value 
of r makes the potential deeper and thus enhances the 
growth of non-axisymmetric modes, with a slightly larger 



growth rate, F ~ 1.41t d 



5.9 Gyr for r = id/24 versus 



i.33t: 



5.6 Gyr -1 for r -» 0. 



Figure [Tp shows the ratio 



o(°) 



s(2) 



B 



(0) 



(42) 



where the last equality does not hold for numerical solutions, 
which include higher-order modes, and where, 



r(m) 



B<f, 
Da, 



(43) 



The quantity 5 (2) is given by the second term in the brack- 
ets of Eq. ((29} . The figures show 5 at the azimuth (j> = 
0cor — Kjfcor /2 where one of the a-spiral arms crosses the co- 
rotation circle. The (^-independent) amplitude (envelope) of 
this ratio is shown with dashed curves, and a(r, <^w) is rep- 
resented with a dash-dotted curve for reference. For both 
finite and vanishing r, the envelope of 8 exceeds unity near 
the co-rotation, so that the jmj = 2 modes dominate there. 
The dominance is stronger when r differs from zero. Impor- 
tantly, the radius r max where the envelope of S is maximum 
is somewhat larger for the finite r. 

The azimuthal phase differences A^ and A r between 
the maximum in, respectively, the azimuthal and radial mag- 
netic field components, and the maximum of a, are shown 
in Fig. QJ. For both r — > and t / 0, the magnetic spi- 
ral arms (whose phase shift is given approximately by A^, 
since B$ dominates over B r ) cross the a-spiral near the co- 
rotation radius, so that each magnetic arm leads the a-arm 
in azimuth inside the co-rotation circle and lags it at larger 
radii. This implies that the magnetic arms are more tightly 
wound than the material arms. 

The r effect produces a phase shift in the magnetic arm 
such that the part of the arm that leads the corresponding 
a-arm is weakened, while the part that lags the a-arm is 
enhanced. This can be seen from Fig. For the r — > 
case, r max (vertical dotted red line) is located slightly inside 
the circles at which = and A r = 0. At these circles, 
the Ba, and B r magnetic arms respectively cross the a-arm. 
This means that for r — > 0, the part of the magnetic arm 



that leads the a-arm is stronger than th at which lags. (For 
the more realistic numerical solutions of |Paper~il . the lead- 
ing and lagging parts are of equal strength for t -> as 
can be seen in Fig. 8c of that paper.) On the other hand, 
for r = 21 /u in our asymptotic solution, we see that r max 
is located to the right of the radii for which A^ = and 
A r = 0, which tells us that the lagging part of the magnetic 
arm is stronge r than th e leading part. This is consistent with 
the findings o f IPaper J . Thus, the lagging, outer part of the 
magnetic arm is enhanced for finite r, while the leading, in- 
ner part is weakened. For example, the phase shift between 
the magnetic and a-arms is A^(r ma x) = 5.2° for r — > but 
—3.7° for r = 21 1 u. However, what we are most interested 
in is the effect of a finite r, so from now on we will mostly 
consider the phase difference between the vanishing and fi- 
nite r cases. Crucially, we shall show below in Sect. I4.3l that 
the general behaviour of the phase difference as a function 
of fir is nicely reproduced by the asymptotic solution. Thus, 
defining the difference between the phase shifts obtained for 
r = and finite r, 



A(r) = A^ [r, r max (r)] - A [0, r max (0)] . 



(44) 



we have A(2Z/u) = —9° for the asymptotic solution, which is 
of the same order of magnitude as —fir = —14°. The phase 
differenc e A(r) ~ —Sir, obtained here and in the models of 
IPaper j . is a natural consequence of the finite response time 
r of £ to variations in a. 

In Fig. [JJi, we plot the pitch angle of the magnetic field 
Pb at the same azimuth (j> = <f> COT — Kr COI /2, shown with 
the opposite sign for presentational convenience. The pitch 
angle obtained for the axisymmetric disc (dashed lines) is 
about pb — —(20-30)° in the region shown, and increases 
with radius in agreement with Eq. (|37[) . Its large magnitude 
is due to the large, constant value of ao used in this model, 
whereas the increase of |j>b| with r is due to our choice of 
a constant disc scale height here, as discussed above. The 
azimuthal variation in a leads to a large variation in the 
pitch angle, with — ps being large where a is large and small 
where a is small (compare with the dash-dotted curve for a 
in Fig. [Tp). The magnitude of the variation depends on e a . 
This variation is explained by the fact that the auj dynamo 
mechanism tends to produce \B r /B$\ < 1 due to the veloc- 
ity shear; hence, \ps\ oc \ao/G\ ' 2 , Eq. (f35)l . (By contrast, 
the a 2 -dynamo would generate rather open magnetic spirals 
with \B~r/B~4\ ~ 1.) 



4 COMPARISON WITH NUMERICAL 
RESULTS 

4.1 Numerical solutions of Paper I 

The semi-analytical solution obtained agrees w ell, qualita- 
tively, with the numerical solutions of IPaper J. whi ch can 
be seen by comparing Fig. [Q with Fig. 8 of IPaper 1 which 
shows the solution of the fiducial non-axisymmetric model of 
that paper (Model E) in the kinematic regime. This is espe- 
cially true in the r — > case. The numerical solution, how- 
ever, does reveal a greater difference between the finite and 
vanishing r cases than what is seen in the semi-analytical 
solution . For example, the phase difference for Model E of 
IPaper J during the kinematic stage is A(//«) = —27° ± 4° 
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Figure 2. The relative strength S of the non-axisy mmetric part 
of the mean magnetic field in Model E of lPaper~3 . defined as in 
Eq. (142 1 . as a function of the galactocentric radius near the co- 
rotation circle at <j> = <f> C or (black solid), and the envelope of this 
dependence at all values of <j> (black dotted). The individual con- 
tributions of the \m\ = 2 [red dashed for <5' 2 ) (r cor ) and red dotted 
for its envelope] and \m\ = 4 [blue dash-dotted for <5' 4 '(r C or) and 
blue dotted for its envelope] modes are also shown, (a) r — > 0; 
(b) r = l/u. 



where the range of A reflects the finite resolution of the nu- 
merical grid. This is comparable to, though somewhat larger 
in magnitude, than —Sir = —18° for that model, while in 
our semi-analytical solution, presented above, A was found 
to be comparable, but somewhat smaller in magnitude, than 
—Q.T for the disc model used. This apparent discrepancy sug- 
gests that there is an additional physical element missing in 
the asymptotic solution. In fact, the larger phase shift in the 
numerical solution can be partly attributed to higher-order 



even corotating modes, especially m = ±4. Figure [2^, shows 
the relative strength of the non-axisymmetric magnetic field 
8, defined in Eq. (|42|l . at the azimuthal angle <f> = C or, where 
the a-spiral cr osses the co-rotation circle, for the case r — > 
in Model E of |Paper"H Fig. [2p shows <5(0 C or ) for the case 
r = l/u. These panels are similar to Fig. 8b of lPaper J . with 
the envelope of 8 over all cj> plotted here as a dotted line. 
The quantities <r 2 ' and <r 4 ' [see Eq. [43] , are also plotted. 
It is clear that the \m\ — 2 modes dominate, as assumed in 
the asymptotic solution. However, it is also clear that the 
ratio of the amplitude of \m\ = 4 to that of \m\ = 2 is larger 
when r is finite. Moreover, for r = l/u, the <5' 4 ' envelope 
peaks at a larger radius than that of <r 2 ' . This suggests that 
the \m\ — 4 modes have the effect of shifting r max outward. 
Since A(r) is sensitive to r max (r), this can cause a significant 
increase in A(r). 

Although one could, in principle, include modes with 
7/1 1 > 2 in the asymptotic solution presented above, every 
extra mode would add four new differential equations as well 
as new terms to the existing equations (|16|l - (|19[) . greatly 
increasing the complexity of the model. We take the view 
that the asymptotic solution presented above does a very 
good job of reproducing the key qualitative features of the 
numerical solution, and that further refinements for the sake 
of numerical accuracy are not necessary at this time. 

4.2 Direct comparison of asymptotic and 
numerical solutions 

It is also useful to compare an asymptotic solution such as 
that presented above with a numerical solution for the same 
disc model ( with e.g. a much smaller half-thickness h than 
the models of lPaper J ). An example is shown in Figs. [3] and 
13] for t — ¥ and r = l/u, respectively. The model uses the 
same parameters as the asymptotic solution above, except 
that e = 0.6 instead of 1, for both the asymptotic and numer- 
ical solutions. Plots are shown at t = 10id, when the relative 
strength of higher-order radial modes has become quite low 
(though they can still be seen as wiggles at r > 11 kpc). 
In Figs. [3^, (asymptotic) and c (numerical), we compare the 
quantity 8 for the two types of solution, while in[3p (asymp- 
totic) and d (numerical) , the quantities A r and A^ are com- 
pared, for the t — > case. For this case of vanishing r, the 
two types of solution are in fairly good agreement. However, 
as can be seen by comparing panels a and c, magnetic arms 
are more pronounced in the numerical solution than in the 
asymptotic solution. This can be attributed to higher order 
enslaved modes in the numerical solution; as seen by com- 
paring the overall envelope (dashed) with the envelopes of 
the |m| = 2 (dotted) and |m| = 4 (dash-dotted) modes in 
panel c. 

For the t = l/u case, shown in Fig. [4] qualitative fea- 
tures of the solution are in reasonable agreement, but the 
effect of a finite r is more dramatic in the numerical solu- 
tion than in the asymptotic solution. For instance, r max is 
much larger in the numerical solution, resulting in a larger 
phase shift. This is seen by comparing panels b and d, where 
it is evident that the black dotted vertical line crosses the 
black solid and dashed lines at much larger values of |A^| 
and | A r | in t he numerical solution. As with the solutions 
from lPaper J . enslaved modes with \m\ > 2 are stronger in 
the finite r case than in the r — > case, as can be seen 



© 2012 RAS, MNRAS 000,[THT4] 



10 L. Chamandy, K. Subramanian and A. Shukurov 




Figure 3. Direct comparison of asymptotic (a-b) and numerical 
(c-d) solutions for the same disc model, for the case r — > 0. Panels 
a and c are similar to panel b of Fig. [TJ whilst panels b and d are 
similar to panel c of that figure. In a and c the \m\ = 2 and 
\m\ =4 envelopes are represented by dotted and dash-dotted 
curves, respectively. In Panel d, the curves are smoothed over 
0.125kpc. 

by comparing the dash-dotted envelopes for the \m\ = 4 
modes in Figs. 0fc and (3^. The assertion, made above, that 
higher-order enslaved modes are responsible for enhancing 
the phase shift is therefore strengthened. 

The fact that non-axisymmetric modes can be easier 
generated in a thinner disc, where the difference in the global 
mode structure is less important because of the stronger 
dominance of the local dynamo action (due to the shorter 
magnetic diffu sion tim e across the disc), is well known 
(Sect. VII. 8 in RSS88). However, the importance of the fi- 
nite relaxation time in this respect was not appreciated ear- 
lier. We emphasize that this disc model is used for the sake 
of illust ration only, and is not as realistic as the models of 
iPaper J . where enslaved modes with \m\ > 2 are anyway 
much less important. 

4.3 The phase difference caused by a finite 
dynamo relaxation time 

It is also interesting to explore how the phase difference A 
varies with the dimensionless quantity Sir. Here we retain 
the d isc mode l (of the asymptotic example or of Model E 
from IPaper Jl but vary r. Figure [5] shows the phase dif- 
ference A as a function of Sir for the asymptotic solution 
of Sect. 13.31 shown by crosses in the plot, confirming that 
A ~ —CQt with C a constant of order unity. It can be seen 
that the relation is not strictly linear and tends to flatten as 
Sir increases, but is approximately linear for Sir < 0.5, which 
is the region of parameter space that normally applies to 
disc galaxies. Figure [S] also pr esents re sults from the numer- 
ical solution of Model E from lPaper II . discussed in Sect. 14. ll 
above, where the galaxy model is more realistic than in the 



Figure 4. Same as Fig. \3\ (but note the change in the plotting 
range of <5), now for r = l/u. 



analytical solution obtained in Sect. 13.31 For Model E, the 
scaling relation A oc Sir (for small Sir) quickly establishes 
itself in the kinematic regime, as soon as the leading eigen- 
function becomes dominant. Remarkably, it applies to both 
kinematic dynamo solutions and to the steady state, with 
more or less the same proportionality constant. The magni- 
tude of the proportionality constant C « 1.7 (represented 
by the solid line in the figure) is still of order unity, but 
is significantly larger than C ~ 0.65 (dashed line) found in 
the asymptotic solution. Although this difference may be 
partly attributable to the different disc models used, the 
larger value of C in the numerical model is also partly due 
to higher-order even modes, as discussed in Sect. 14. ll above. 

In summary, asym ptotic so lutions obtained here clar- 
ify numerical results of |Paper~H and help us to isolate their 
generic features. Figure [T] can be compared with Fig. 8 of 
IPaper J . which shows similar kinematic results, confirming 
the presence of a strong \m\ = 2 magnetic component near 
r — r COI . Around this radius we see clearly from the asymp- 
totic analysis that the effective potential V(r) has a new 
minimum, allowing new bound states (or growing modes) of 
q(r) near co-rotation. Both asymptotic and numerical solu- 
tions reveal stationary even modes in the frame corotating 
with the spiral pattern of the a-coefficient. Both have mag- 
netic spiral arms more tightly wound than the a-arms, with 
the two types of arm crossing near the co-rotation circle. 
Moreover, both numerical and asymptotic results, the latter, 
e.g., in Eq. (J3U), also agree in that the response of the mean 
magnetic field to the a-spiral pattern is stronger when r is 
finite. More importantly, they agree in that magnetic arms 
undergo a negative phase shift of order —Sir such that the 
part of the magnetic arm which trails the a-arm is enhanced 
when r is finite. Enslaved modes with \m\ > 2, though weak, 
can substantially enhance the phase shift because they are 
located further out in radius than the \m\ = 2 mode, where 
magnetic arms lag a-arms. 
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Figure 5. The phase difference A between the solutions with r ^ 
and t — as a function of Qt. Results for the parameter values 
used in the illust rative exa mple are shown with blue crosses; those 
from Model E of iPaper"!! are shown during the kinematic regime 
(at t = 3td o after the simulation has begun) and in the steady 
state. 



5 CONCLUSIONS 

We have extended non-axisymmetric mean-field dynamo 
theory in two important ways: 

• Obtained semi-analytically a solution for the axisym- 
metric and strongest enslaved non-axisymmetric modes 
(m = and m = ±2 for a two-armed spiral). Previously, 
the only analytical treatment that existed was for the sub- 
dominant non-enslaved modes (m — ±1); 

• Included the effects of the finite relaxation time of the 
mean electromotive force (the r effect) in our solutions. 

The asymptotic solutions employ the WKBJ approxi- 
mation applicable to tightly wound magnetic spirals; to ob- 
tain explicit expressions for the growth rate of the mean 
magnetic field and its radial and azimuthal distribution we 
use the no-z approximation (known to be quite accurate for 
thin discs) to solve the local dynamo equations. 

We find over all good qualitative agreement with the nu- 
merical results o f lPaper J for a global, rigidly rotating spiral 
pattern (which is the type of spiral forcing most amenable to 
analytical treatment) to grow around the co-rotation radius. 
In particular, the asymptotic s olution o btained here agrees 
with the numerical solution of lPaper II in the following key 
features. 

(i) Strong magnetic modes with the azimuthal wave num- 
bers m = ±n are supported by the dynamo action, which 
corotate with the a-spiral (here n is the number of a-spiral 
arms). Non-axisymmetric modes are confined to a radial 
range of a few kiloparsecs around the co-rotation radius and 
are comparable in strength with the m — mode. 

(ii) Magnetic arms (understood as ridges of the mean 
magnetic field strength that have a spiral shape) cross the 
a-spiral near the co-rotation radius and are more tightly 
wound than the a-arms. 

(iii) In the model considered here, the only effect of the 
galactic spiral pattern is the enhancement of the dynamo 
Q-coefficient by the spiral pattern. As a result, the mag- 
netic lines of the mean field are less tightly wound (larger 
magnitude of the pitch angle pb) within the a-spirals than 
between them. 



(iv) The amplitude of the non-axisymmetric modes in- 
creases with the magnitude of the relaxation time r. 

(v) Another effect of the finite magnitude of r is that the 
part of the magnetic arm that lags behind the a-arm with 
respect to the galactic rotation (i.e., outside the co-rotation 
radius) becomes stronger than the part that leads the a- 
arm (inside the co-rotation radius) . This leads to the overall 
appearance of magnetic arms that are phase-shifted from 
the a-arms (material arms) in the direction opposite to the 
galactic rotation. 

(vi) For fir < 0.5, the phase shift produced by a finite 
value of r (relative to that for r — > 0), is directly pro- 
portional to —Qt, with a proportionality constant of order 
unity. 



6 DISCUSSION FOR PAPERS I AND II 



Results obtained in IChamandv et all (|2013T ) dPaper j) and 
here (Paper II) for a rigidly rotating a-spiral demonstrate 
the ability of the mean-field dynamo mechanism to produce 
magnetic arms that do not overlap with the material arms 
in the regions where the former are best pronounced, which 
is outside the co-rotation circle in our model. Due to a finite 
magnitude of the relaxation time of the mean electromotive 
force, t , the lagging part of the magnetic arm is better pro- 
nounced than that ahead of the material arm, which can 
make it difficult to detect their intersection in the observa- 
tions. As a result of advection of the mean magnetic field by 
the differentially rotating gas, magnetic arms are, in these 
models, more tightly wound than the a-spiral which drives 
them, and are localized to within a few kiloparsecs of the co- 
rotation circle. The phase shift between the magnetic and 
a-arms in our model vanishes at or near the co-rotation ra- 
dius and then increases, first rapidly and then at a lower 
rate, up to a value as large as 60° at a distance of a few kpc 
from the co-rotation circle. 

The morphology of magnetic arms in spiral galaxies 
and their position relative to the material arms is rather 
diverse. One extreme is the galaxy NGC 6946 where mag- 
netic arms appear to have the same pitch angle as the optical 
arms, and the two patterns do not intersect. The azimuthal 
phase shifts between the individual optical arm s and their 
magnetic counterparts identified by iFrick et all (|200Ch arc 
-26° ± 12°, -36° ± 11° and -45° ± 17° at r « 4 kpc and, 
for an outer arm, —68° ± 17° at r as 8 kpc, where negative 
phase shift corresponds to a position behind a gaseous arm 
(all distances ha ve been reduced to a distance of 5.5 Mpc 
to NGC 6946 - iKennicutt et al.ll2003l). Accordin g to the 
Ha observations and analysis of Fathi et alj l|2007h . the co- 
rotation radius of the outer spiral pattern in this galaxy is 
at about 8.3 kpc, whereas the interlaced magnetic and opti- 
cal spiral arms occur at 1 < r < 9 kpc. The dynamo action 
at r > 9 kpc may be too weak to make the magnetic arms 
observable there, but there are no indications of the mag- 
netic and optical arms crossing at or near the co-rotation 
radius. Nevertheless, the right magnitude of the azimuthal 
phase shifts obtained in our models (15-40°) is arguably 
encouraging. 

Another possibility is that the spiral structure in 
NGC 6946 is more complex than that of a single, rigidly 
rotating pattern. For instance, the arms may be winding 
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up t o some extent, and thus transient fe.g. iDobbs et al 



20ld; ISellwoodll201ll ; lOuillen et al.ll201ll; IWada et al.ll201ll; 
Grand et all 120121 ). The model of IComparetta fc Quiller] 



1 20121 ) invokes multiple rigidly rotating patterns, which in- 
terfere to produce spiral features that wind up. Magnetic 
arms which are present over several kiloparsecs in radius, 
and which have a large negative phase shift from the a-arms 
that varies only weakly with radius, like in NG C 6946, a re 
indeed found in our 'winding-up' spiral model of iPaper 1 

On the other hand, the mutual arrangement of mag- 
netic and material arms is more comp licated in M51 
[|Patrikeev et al.ll2006l ; iFletcher et al.ll201ll ) where the two 
patterns overlap in some regions but are systematically off- 
set by about 0.5-0.6 kpc elsewhere (adopti ng 7.6 Mpc for the 
distance to M51 - Idardullo et al.ll2002l ). This linear dis- 
placement corresponds to the angular phase difference of 5- 
10° at r = 3-6 kpc where such an offset is best pronounced. 
The two arms visible in polarized emission at A6.2cm inter- 
sect t he material arms obser ved in other tracers at r = 5- 
6 kpc. Elmegrc en et al.l |l989) suggest that M51 has two spi- 
ral arm systems, with the co-rotation radii at r = 4.8 kpc 
and r = 12 kpc. iGarcfa-Burillo et al.l i|l993l ) determine the 
co-rotation radius to be 5.8 kpc from CO observations and 
modelling of M51. Then the region at r = 3-6 kpc is (mostly) 
inside the co-rotation radius and the magnetic arms in the 
inner galaxy are displaced downstream of the material arms, 
as expected in our model. The two magnetic arms are posi- 
tioned differently with respect to the material arms at larger 
galactocentric distances. Outside the co-rotation circle, mag- 
netic Arm 1 (on the east of the galactic centre) is system- 
atically lagging the material arm, in accordance with our 
model, but the magnetic Arm 2 overlaps the material arm. 

Our model appears to be better applicable to this 
galaxy, but the diverse mutual arrangement of the mag- 
netic and material arms in M51 strongly suggests that 
more than one ph ysical effect is involved in its genesis. 
IDobbs et al.l (|2010h suggest that the spiral arm morphol- 
ogy in M51 evolves rapidly due to the interaction with its 
satellite galaxy. Mag netic arm s driven by evolving material 
arms are explored in |Paper~H It also appears that M51 may 
be a good candidate for the 'ghost' magnetic arms, where 
the magnetic arms trace the spiral arms as they were in 
th e past up to a few h undred Myr earlier. The simulations 
of IDobbs etall (|2010h show that the material spiral arms 
are being wound up by the galactic differe ntial rota tion, so 
the model with transient material spirals of iPaperll may be 
more relevant to M51 than the model with a rigidly rotating, 
stationary spiral pattern explored here. 

Still another morphology of magnetic and gaseous 
arms is observed in the nearby barred galaxy M83. 
iBeck et al.l (|in preparation!) ap plied the same ap proach in- 
volving wavelet transforms as [Prick et all (|2000l) to deter- 
mine the positions of the spiral arms as seen in optical light, 
dust, Ha, CO, as well as the total and polarized intensi- 
ties at AA6 and 13 cm. The morphology of M83 is domi- 
nated by a well-pronounced two-armed spiral pattern in each 
tracer, with some substructure within the arms. The mag- 
netic and material arms in this galaxy are clearly separated 
but intersect at the galactocent ric radius 6.4 kpc (given the 
distance to M83 is 4.5 Mpc - iThim et al.ll2003h. The co - 
rotation radius of M83, determined by Hirota et alj (|2009l ) 
as 2.4 arcmin, or 6.2 kpc, is practically equal to the intersec- 



tion radius given the uncertainties involved. This galaxy is 
a good candidate for the formation of magnetic arms by the 
mechanism suggested here. 

An important assumption of all or most of the avail- 
able models of magnetic arms is that the a-effect is stronger 
within, or at least correlated with, the material spiral 
arms. We note that the effects of the spiral arms on 
the scale height of the interstellar gas, its turbulent scale 
and velocity, local velocity shear, and other parameters 
that control the intensity of the dynamo action are far 
from being certain, either observationally or theoretically 
(see IShukurovlll998l; IShukurov et al l |2004| . and references 
therein). IShukurov fc Sokolofj 1 19981) argue that the mag- 
nitude of the a effect within the material arms can be four 
times smaller than between the arms, whereas the turbulent 
diffusivity is plausible to be only weakly modulated by the 
spiral pattern. As a result, the local dynamo number within 
the material arms can be four times smaller than between 
them, and the dynamo action can, in fact, be more vigorous 
between the material arm s. Direct numerical simulations by 
lElstner fc Gressell l|2012l) find that both the a effect and tur- 
bulent diffusivity increase with star formation rate a, which 
itself will be larger within the arms, but this leads to an 
overall dynamo number that scales inversely with a. How- 
ever, even if the dynamo number is larger in the interarm 
regions, this does not necessarily imply that the steady-state 
mean magnetic field there should be stronger than within 
the arms, since the gas density is lower in between the arms 
than within them. Significant effort in observations, theory 
of interstellar magnetohydrodynamics and dynamo theory 
are still required to establish a clear understanding of the 
various mechanisms that produce the diverse morphology 
of magn etic arms in spiral galaxies. Our work here and in 
iPaper il provides the beginnings of such a study. 



APPENDIX A: THE LOCAL 
EIGENFUNCTIONS 



With approximations described in Section 13.2.11 Eqs. (|18[1 
and (|19p yield, after some algebra, 



b ±2 = b ^\D \X ±2 , a ±2 = -Y ±2 b ±2 . (Al) 
ao 



Here 



and 



A 2 + B 2 ^A 2 + B 2 



= 2, 



(A2) 



Y ±2 (r) = ATiB=\/A 2 + B 2 e^=Y; 2 , (A3) 
with the dynamo number defined as 



a Gh 3 aoGt\ 

L>0 = o = 1 < U ' 

V? h 



(A4) 



where t& = h 2 /rjt is the vertical turbulent diffusion time 
scale and 



R A ■ R B 

cos p — — , smp : 



VA 2 + B 2 



VA 2 + B 2 



(A5) 
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with 

+ [(i + rv) 2 + (2ttT) 2 ] 1 

x [4(i + rv)(±7r 2 + ^V)^ 2 + e' b 2 h 2 ) 

- 2c 2 T {ir + 8' a 2 h 2 ) e'X + (in 2 + e' 2 h 2 ) e'X] 

-4(i + rr)0^ 41 



+ cvrta [W + e' 2 h 2 + e' b 2 h 2 ) ~ 2c T zt d h 2 (6" + > 



B = —I - 2firt 2 (r 2 - 4£ 2 ) + 4(1 + rr)rwid 
2c T I 

+ [(l + rr) 2 + (2fir) 2 ] _1 

X [2c 2 T Qr{\n 2 +6' 2 h 2 ){\-K 2 + e' 2 h 2 ) 



(A6) 



+ c 2 (i + iv) [(Itt 2 + e ,2 h 2 )e' b 'h 2 + {\n 2 + e ,2 h 2 )e':h 2 ] 

- 2c 2 T Q.re"e b h 4 ^ 

+ 2c T st d (iyr 2 + e' 2 h 2 + e' b 2 h 2 ) + orTt d h 2 (e'a + o't 



Further, 

cos /3 

with 
A 



A 



\J A 2 + B 2 



, sin/3 



B 



\J A 2 + B 2 



r (I + Ft^tt 2 + e' b 2 h 2 )-2QT9' b 'h 2 



G [(1 + Yt) 2 + (2nr) 2 ]Gt A 

2w , _ 2SIt{\-k 2 + 6' b 2 h 2 ) + (l + TT)6' b l h 2 



(A7) 



(A8) 



(A9) 



(A10) 



[(l + rr) 2 + (2fir) 2 ]Gi d 
Substituting fl5)l. (|A2|> and (|A3)l into Eq. (|Atj) . wc arrive at 



a±2 = — - 



b±2 



D \b e a A 2 + B 2 



2 V A 2 + B 2 
I D 1 boe a 1 



2 V^4 2 + B 



exp [=Fi(«r + P + P)\ , (All) 
= exp[=R(«x + /3)]. (A12) 
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